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Abstract. Data for the broadening of 24188 Fe II lines by collisions with neutral hydrogen atoms have been computed 
using the theory of Anstee & O'Mara as extended to singly ionised species and higher orbital angular momentum states by 
Barklem & O'Mara. Data have been computed for all Fe II lines between observed energy levels in the line lists of Kurucz with 
log gf > —5 for which the theory is applicable. The variable energy debt parameter E p used in computing the second order 
perturbation theory potential is chosen to be consistent with the long range dispersion interaction constant Ce computed using 
the /-values from Kurucz. 
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1. Introduction 



The line broadening theory of Anstee & O'Mara ( fT99T1 [T995t 
was developed for s-p and p-s transitions between low-lying 
states of a neutral atom by neutral hydrogen collisions, with 
the hope of being able to compute data for a large number of 
lines of various elements for astrophysical applications. The 
approximations necessary for development of such a gener- 
ally applicable theory led to an expectation that the theory 
might enable cross-sections to be calculated with accuracies 
of around 20%; a significant improvement on the classical van 
der Waals description (Unsold 1955 1, commonly used in stellar 
spectroscopy due to its general applicability, with uncertain- 
ties of typically a factor of 2 or 3. Application of the data to 
solar lines of a number of elements, including Na, Ca and Fe 
(Anstee & Q'Mara [T995l Anstee et al. 119971 . in fact suggested 
that the line broadening cross-sections might have uncertain- 
ties as low as 5% for the considered lines, as judged from the 
internal consistency and agreement of abundances with the me- 
teoritic abundances. Anstee & O'Mara ( 1991 1995 ) suggested 
that this may be the result of error cancellation. Recent com- 
parisons with detailed calculations (eg. Kerkeni et al. [2004 
Barklem & O'Mara 2001) indicate uncertainties in the range 
5-20%, the largest differences occuring in lines where the in- 
fluence of avoided ionic crossings turns out to be particularly 
important in the adiabatic interaction potential, an effect which 
is neglected in the Anstee & O'Mara theory to retain generalil- 
ity. The theory has been extended to transitions between higher 
orbital angular momentum states (Barklem & O'Mara 1997 
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Barklem et al. 11998b and to the case of a singly ionized atom 
(Barklem & O'Mara [T9981 . We will refer to these collected 
works by Anstee, Barklem and O'Mara as ABO theory. The 
extended theories have also been tested against the solar spec- 
trum where possible and there is no evidence that these calcula- 
tions have uncertainties significantly different from the earlier 
calculations. 

The theory for ionized atoms differs significantly from the 
theory for neutral atoms, in that calculations must be done line 
by line. In the ABO theory, where potentials are computed us- 
ing Rayleigh-Schrodinger perturbation theory, a key parameter 



is E p , the fixed energy debt in the Unsold ( 1927 1 approximation 
to the second order interaction energy between the perturbed 
absorbing atom and the ground state H atom. For neutral atoms 
energy level splittings are typically small compared to those 
for the hydrogen ground state and thus hydrogenic contribu- 
tions are expected dominate the long range interaction. In this 
approximation, also due to Unsold ( 19551, E p « —4/9 atomic 
units irrespective of the species or state of the neutral atom. 
This approximation enabled general tables of line broadening 
data to be computed, which could be interpolated in for any 
species. In ionized atoms, energy spacings are typically larger 
than those in neutral atoms, and thus the Unsold approxima- 



tion of E p = -4/9 is questionable. Barklem & O'Mara ( 1998 1 
showed that while for the considered lines of Ca II and Ba II 
this approximation was reasonable, for the resonance lines of 
Mg II it led to an error of over 30%. Such a large errror, when 
combined with the intrinsic error in the theory, is considered 
unacceptably large for astrophysical application. Thus, for ion- 
ized atoms calculations would need to be done on a line to line 
basis. 

Up to now, calculations for ionized atoms have only been 
performed for a small number of astrophysically important 
lines (Barklem & O'Mara 1T998I I2000i due to the need to 
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choose a suitable value for E p . The value of E p may be in- 
ferred from the dispersion (i.e. van der Waals) component of 
the long range interaction Cq/R 6 , by choosing E p such that 
the potential will have this long range behaviour, and assum- 
ing this value at all separations. Calculation of C@ involves a 
summation over all /-values connected to the states of the two 
atoms. In these past works the /-values for the perturbed atom 
have being collected from the literature, making the procedure 
labourious. For some species, an alternative approach to this 
problem is provided by the line lists computed by Kurucz us- 
ing a semi-empirical least-squares approach to Slater-Condon 
theory (Kurucz [T973l Cowan [T98"Tl ?16-3). While the /-values 
computed via this method are perhaps not as accurate as ex- 
perimental values which might be found in the literature for a 
small number of selected lines, the enormous advantage is an 
unparalleled degree of completeness which is important for our 
problem. As we will discuss, the accuracy should be sufficient 
for our purposes. The line lists provide not only the required 
/-values for computing Cq, but, of course, an extensive list of 
lines for which we can compute line broadening cross-sections. 
We have computed such data for Fe II, and the calculations and 
results are now described. 

2. Line Broadening Calculations 

The theory, particularly the potentials and collision dynamics, 
has been described in the earlier ABO theory papers. For our 
specific case of a singly ionised Fe atom, the theory is described 
particularly by Barklem & O'Mara J 19981 hereafter B098, and 
references therein). 

To perform the calculations we have used information 



from three files from Kurucz (2003; kurucz.harvard.edu): the 
Fe II energy level list (gf2601.gam), the complete line list 
(gf2601.1in), and the line list (gf2601.pos) which includes only 
lines between observed levels. These files were computed in 
2003, following the method earlier described by Kurucz ( 119731 
1981 1 but with updated observed energy levels and identifica- 
tions. 

First, the various input parameters for the ABO theory must 
be adopted for each energy level in the Kurucz energy level list 
(gf2601 .gam). The ABO model for the interaction potential de- 
scribes the absorbing atom as a positively charged core with a 
single optical electron. The unperturbed optical electron is con- 
sidered as moving independently of the core electrons in a cen- 
tral field, and thus its state may be characterised by principal 
quantum number n and orbital angular momentum quantum I 
along with its binding energy. Together, these three pieces of 
information allow the unperturbed wavefunction for the optical 
electron to be computed in this approximation. Such a descrip- 
tion is justified by the fact that the broadening is dominated by 
separations on the tail of the electronic wavefunction (Anstee 
& 0'Mara [T99lt . 

For each level we determine this information from the dom- 
inant electron configuration and term provided in the Kurucz 
energy level list. In all cases we assume the electron (or one 
of the electrons) with the highest I of those electrons with the 
highest n in the dominant configuration is the optical electron, 
since this electron will have the most extended wavefunction. 



The binding energy of the electron is given by the difference 
between the energy of the state E and the series limit i?ii m ;t 
which defines the Fe III parent state. The idea of parentage is 
only clear in the case of, like our model, a single running elec- 
tron outside a more tightly bound core. In more complex sit- 
uations, particularly where there are multiple open shells, one 
often has parental mixing (see eg. Condon & Odabasi 1980 
Cowan frWTV 

The most common electron configuration in Fe II is of 
the form 3d 6 (a' Mp L p )nl a M L where a' Mp L p defines clearly 
the parent term in Fe III. However, in other configurations the 
parentage may be unclear or parental mixing occurs. Levels of 
the form 3d 5 (a' Ms L g )AsAp a M L, where a' Ms L g defines the 
grandparent, assuming the Ap electron to be our optical elec- 
tron, has the possible parents Mg±1 L g . If M = M g + 2 then 
we must have M9+1 L g , while if M — M g — 2 we must have 
Mg ~ 1 L g . If M = M g then either possibility may occur and 
we arbitrarily chose the alternative giving the lowest binding 
energy, and thus the most extended wavefunction. 

While our model is not really applicable to the case where 
the external shell contains equivalent electrons, we proceed 
on the assumption that a reasonable estimate might be ob- 
tained by taking one of the equivalent electrons as the opti- 
cal electron. We note that since the Cq calculation is not de- 
pendent on the optical electron model, our long range poten- 
tial will be correct since it is referenced to the Cq calculation. 
For the configuration 3d 7 , which has a large number of frac- 
tional parents, we used the principal parent (see Condon & 
Odabasi 1980 1, namely a 3 F. The configuration 3d 5 4s 2 a M L 
has two possible parents, namely M±1 L. In an attempt to best 
model the wavefunction of the 4s electron with the most ex- 
tended wavefunction we assigned the parent as that parent 
which gives the smallest binding energy which is M+1 L by 
Hund's rule. Finally, for the small number of levels of the form 
3d 5 {a' A ^L g )Ap 2 ( 1 S, 3 P, 1 D) a M L treated in this work we 
have assumed that the parent is z"' P unless M ~ 4 in which 
case we assume z 5 P. 

Having assigned to each state n, I and parent terms, and 
thus binding energies for the optical electron, scaled Thomas- 
Fermi-Dirac radial wavefunctions (Warner 1968 1 were then 
computed for the optical electron for each state. Due to neglect 
of exchange effects ABO calculations are only valid for low- 
lying states where such effects are not large (see relevant ABO 
papers), and thus we treat only states for which the effective 
principal quantum number n* is less than 3 for s and p states 
and for which n* < 4 for d states. 

Now we need to specify E p for each state. For a perturbed 
atom A with states k' in given energy level k, interacting with 
a H atom in the Is ground state, denoted by I, with excited 
states the Cq value is computed via the standard expression, 
in atomic units, 



EE 



fA rH 
Jkk'JW 



2 t^k V^l ( AE k'k 



(1) 



J k'k L 



where AE*, m — E*, —E*. Data for bound-bound transitions 
in the Fe II atom are taken from the Kurucz complete line list 
(gf2601.1in). The bound-bound transition data for H are taken 
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from Wiese, Smith and Glennon ( 1966|l. As in B098, bound- 
free transitions are accounted for in the H atom, but not in the 
perturbed atom. In B098 an approximate constant Gaunt factor 
for the bound-free component was adopted; a polynomial fit to 
the calculations of Karzas & Latter ( 1961 1 is now used. Once 
the Cq coefficient is computed, the value of E p may be chosen 
such that the ABO theory Rayleigh-Schrodinger- Unsold poten- 
tial will have the correct long range behaviour via 



2(p1) 
Cq ' 



(2) 



The mean square position of the optical electron is com- 
puted from the wavefunction. 

Kurucz has also computed Cq values, provided in his en- 
ergy level list (gf2601.gam), making the approximation that 
the energy levels in the perturbed atom are much more closely 
spaced than the separation of the H ground state from other H 
states (Kurucz fl981l page 76), which leads to 



<kk' 



k'=£k kk 



(3) 



where amis) is the static dipole polarizability of the H atom in 
the ground state, which is 9/2 atomic units (Dalgarno ll962> . 
As a check on our calculations, we also computed Cq in this 
approximation as well as the /-value sum. The results were in 
good agreement with those computed by Kurucz. Comparison 
of our Cq values from equation ^ with those computed by 
Kurucz indicates that this approximation leads to overestima- 
tion of Cq by a mean factor of 1.23, with a standard deviation 
of 0.47, in Fell. 

The results for Cq and E v are presented in Table^as they 
may be of interest for other applications, for example calcu- 
lation of collisional depolarization cross-sections. The static 
dipole polarizabilities are presented since they are easily com- 
puted and such fundamental data may be of interest, for exam- 
ple, in estimating long range interactions with perturbers other 
than H. Note that we consider only observed energy levels (i.e. 
no predicted levels) and states which meet the criteria for n* 
discussed above for the applicability of the ABO theory. 

Line broadening cross sections were then computed for 
lines in Kurucz's list of Fe II lines (gf2601.pos) with log gf > 
—5, where the upper and lower states are both in Table [0 
Potentials are computed employing the computed wavefunc- 
tions and adopting the computed E p values. We have limited 
the calculations as this mechanism is not important in very 
weak lines and computing times are non-negligible, of the or- 
der of a minute per line on a modern workstation. Note that our 
calculations deal only with transitions where Al = ±1; lines 
not obeying this selection rule are typically weak in any case. 

The computed line broadening data are presented in 
Table |21 We have included sufficient information for lines and 
states involved to be uniquely identified. In table |2j we quote 
the line broadening cross section a for a collision speed v = 
10 4 m s~ 4 and the velocity exponent a as defined by equation 
1 of Anstee & O'Mara ( 1995 1. They also provide an expres- 
sion for calculation of the line half half-width per perturber for 



a Maxwellian velocity distribution from these data, their equa- 
tion 3. In table [2] we quote the line full half- width per unit hy- 
drogen atom density for 10 4 K. The temperature dependence 
of the line width is T~ ~ . We also provide the line broaden- 
ing data computed in the Unsold approximation E p = —4/9, 
namely the cross section 0-4/9 an d velocity parameter 04/9. 
These data are of interest in estimating the sensitivity of the 
data to the value of E p . 



3. Discussion 

We now discuss the results including some statistical proper- 
ties, their impact and reliability. 

First, we discuss in general terms the accuracy and reliabil- 
ity of the calculations. Assuming the completeness of the most 
important data for the Cq sum, the accuracy of the computed 
Cq will depend on the accuracy of the /-values for Fe II, the 
data for H being essentially exact. To estimate the accuracy of 
the Kurucz /-values we have compared them with the exper- 
imental data of Schnabel et al. (2004) in Fig.^ The experi- 
mental data have uncertainties ranging from 0.025 dex for the 
strongest lines, to 0.1 dex for the weakest lines. The compari- 
son has been plotted against / /AE, which reflects the weight 
the transition has in the Cq sum (see equation[2and[3j. Overall 
the data compare well, with no evidence for any systematic dif- 
ferences within the uncertainty in the experimental data. The 
mean difference is —0.023 dex with a standard deviation of 
0.31 dex. However, it is clear that the stronger transitions are 
in the best agreement. As shown in the lower panel of Fig.^ 
these are the most important transitions in contributing to the 
Cq sum, and thus the error in these /-values is more relevant 
for estimating the error in Cq. If we consider only those tran- 
sitions with \og(f / AE) > — 1 the mean difference is +0.036 
dex with a standard deviation of 0. 12 dex. The comparison sug- 
gests that the random error in the /-values for the most impor- 
tant transitions is about 0. 1 dex. Thus, we estimate the error in 
the computed Cq values at around this level, i.e. 26%. A rough 
estimate of the error in the broadening cross section can be 
made via the van der Waals theory, in which the cross section 
is proportional to C 6 . Thus, a 26% uncertainty in Cq leads 
roughly to an uncertainty of 10% in the cross section. 

As seen in Fig. there are a small number of quite dis- 



crepant transitions. For reasons pointed out by Kurucz (119811 
pg 29), as in any such large-scale calculation, data that are 
grossly in error may arise. The expectation is, as seen in the 
Fig. Q that these discrepancies are more likely for weak, and 
thus less important, transitions. We cannot, however, rule out 
such large errors affecting important transitions in our calcula- 
tions. It should also be mentioned that our neglect of bound- 
free transitions in the Fe II atom could be important for some 
higher states. 

Figure [2]plots the distribution of the inferred E p values for 
the considered states. We see that the majority of states have 
E p values close to Unsold's hydrogenic value of —4/9 atomic 
units. Notably, there are a number of states for which a low 
value of \E P \ is inferred. As the second order dispersive inter- 
action is inversely proportional to E p , the interaction energy, 
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Table 1. The computed interaction parameters for the 572 considered levels of Fe II. The complete table is available only 
electronically; a short extract is provided here as a guide to its form. For each level tabulated are: the energy E with respect to 
the ground state, the parity (EVE or ODD), the J quantum number, the principal quantum number n for the optical electron, the 
orbital angular momentum quantum number I for the optical electron, the adopted series limit energy E haut , the mean square 
position of the optical electron the dispersion coefficient C§ and the static dipole polarizability of the level, and the 
energy denominator E p . Note that we use the definition AE = Cq/R 6 , while Kurucz uses Av — /i? 6 ; thus the 

definitions differ by a factor of h, in addition to the difference in units. 



E 


Parity 


J 


n 


I 




(pi) 


c 6 


ad 


E p 


[cm" 1 ] 










[cm" 1 ] 


[a.u.] 


[a.u.] 


[a.u.] 


[a.u.] 


0.000 


EVE 


4.5 


4 





130563.000 


8.53 


29.73 


30.5 


-0.574 


384.790 


EVE 


3.5 


4 





130563.000 


8.57 


29.73 


30.6 


-0.577 



Table 2. The computed broadening data. The entire table is available only electronically; a short extract is provided here as a 
guide to its form. The table presents for each line: the species as described by the packed parameter Z + 0.01 (Z — N) where 
Z is the atomic number and N is the number of electrons, the air wavelength A a ir, the J quantum numbers for the lower and 
upper states J\ ow and J U pp, the energies of the lower and upper states E\ ovr and E upp , the line broadening cross-section a for a 
collision speed of 10 4 m s -1 , the dimensionless velocity parameter a, and the log of the line width (FWHM) per perturber at a 
temperature of 10 4 K. Further, the cross section and velocity parameter computed using the approximation E p = —4/9 a.u. are 
given, CT4/9 and a 4 / 9 . 



Species 


Aair 

[A] 


</low 


Jupp 


-E"low 

[cm- 1 ] 


[cm- 1 ] 


a 
[a.u.] 


a 


log(r 10 4 K /jVH) 
[rad s 1 cm 3 ] 


°"4/9 

[a.u.] 


«4/9 


26.01 


890.122 


4.5 


5.5 


0.000 


112344.110 


326. 


0.319 


-7.634 


163. 


0.200 


26.01 


905.210 


4.5 


5.5 


1872.567 


112344.110 


312. 


0.437 


-7.674 


203. 


0.263 



and thus the line broadening cross section, are particularly sen- 
sitive to the choice of E p value when its absolute value is small 
(see Fig. 1 of Barklem & O , Mara l2000l . Figure plots the 
ratio cr/o^/g against the E p value for the upper state. It demon- 
strates that the breakdown of the E p = —4/9 approximation 
can be important when \E p \ is small. We note that the states 
with small \E p \ often correspond to states where the external 
shell contains equivalent electrons. Small \E p \ values arise in 
such cases due to small values (see equation resulting 
from the compact wavefunction for the chosen optical electron 
in such states. Note that Fig.[5]is somewhat deceptive in depict- 
ing the scatter in this ratio, and thus a corresponding histogram 
of 0-/0-4/9 is given in Fig.|4]which will be discussed below. 



Figure [5] compares our line width results with those com- 
puted by Kurucz using the classical van der Waals broaden- 
ing theory (Uns61d ll955> . Our line widths are typically larger, 
with a mean enhancement factor of 1.37 with a standard devi- 
ation of 0.31. Figure|4]compares our best cross sections a with 
those computed in the Unsold approximation 04/9. As expected 
given that most levels have an E p near this value (see Fig. [2), 
for the vast majority of lines the results are in good agreement. 
However, there are a number of lines where the approximation 
breaks down. As shown in Fig. |5J these cases correspond to 
lines with small \E P \ values. 



4. Concluding remarks 

We have computed line broadening data for 24188 lines of Fe II 
from the Kurucz line lists with log gf > — 5 where the ABO 
theory is applicable. Thus, we have demonstrated the possibil- 
ity for such large-scale calculations using the theory. Similar 
calculations should be done for other species, including neutral 
species. Such data will be of importance in both interpreting in- 
dividual lines in cool star spectra and in opacity calculations for 
modelling cool star atmospheres. The influence of collisional 
broadening due to hydrogen on the total line blanketing and 
line blocking can be important (eg. Short & Hauschildt l2004> . 

We estimate the uncertainty in the computations of the dis- 
persion coefficient Cq as approximately 26% leading to an un- 
certainty of 10% in the cross section (§ 3). Based on our dis- 
cussion in section 1, the intrinsic uncertainty for the theory is 
around 10%. Thus, we estimate that the uncertainty in the cal- 
culations is around 20%. Occasional large errors may occur. 
We plan to do similar calculations for all neutral and singly ion- 
ized species for which Kurucz has computed data. The data will 
be incorporated into the Vienna Atomic Line Database ( VALD, 
Kupka et al. fT999l . 



Acknowledgements. We gratefully acknowledge the support of the 
Swedish Research Council. 



PS. Barklem and J. Aspelund- Johansson: The broadening of Fe II lines by neutral hydrogen collisions 



5 




Fig. 1. The upper plot shows a comparison of /-values 
from Kurucz (2003 1 with experimental data from Schnabel 



et al. (2004 1, plotted against the logarithm of the transition's 
approximate weight in a Cg sum, log(//A_E). The lower plot 
shows histograms of log(//A_E) for the most important tran- 
sitions (upper panel), 10th most important (middle panel) and 
100th most important (lower panel) transitions to each C§ sum 
for all considered states in our calculations. 




Fig. 4. Histogram of the ratio 17/(74/9. 




Fig. 5. Histogram of the line width ratios between this work 
and van der Waals theory as computed by Kurucz. 



Fig. 2. Histogram of E p values. The Unsold approximation 
value of E p = —4/9 atomic units is indicated. 
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